Temperature dependence of current self-oscillations and electric field domains in 

sequential tunneling doped superlattices 
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We examine how the current-voltage characteristics of a doped weakly coupled superlattice 
depends on temperature. The drift velocity of a discrete drift model of sequential tunneling in a 
doped GaAs/AlAs superlattice is calculated as a function of temperature. Numerical simulations 
and theoretical arguments show that increasing temperature favors the appearance of current self- 
oscillations at the expense of static electric field domain formation. Our findings agree with available 
experimental evidence. 
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I. INTRODUCTION 



Manifestations of vertical transport in weakly coupled 
semiconductor doped superlattices (SLs) include elec- 
tric field domain formation, multistability, 
self-sustained current oscillations, and driven and 

undriven chaos. Jl0[ Stationary electric field domains 
appear in voltage biased SLs if the doping is large 
enough. jn| When the carrier density is below a criti- 
cal value, self-sustained oscillations of the current may 
appear. They are due to the dynamics of the domain 
wall separating the electric field domains. This domain 
wall moves through the structure and is periodically re- 
cycled. The frequencies of the corresponding oscillation 
depend on the applied bias and range from the kHz to 
the GHz regime. Self-oscillations persist even at room 
temperature, which makes these devices promising can- 
didates for microwave generation. H Numerical calcu- 
lation of the voltage-doping SL phase diagram shows 
that only static electric field domains are possible for 
high enough SL doping. As the doping decreases, volt- 
age windows where current self-oscillations are possible 
open up. |H| These windows may coalesce into a single 
one as doping is further lowered, and oscillations disap- 
pear below a critical doping value. Since doping is not 
a feasible control parameter, other quantities affecting 
carrier density should be used to observe these behav- 
iors. Feasible control parameters are laser illumination 
in undoped SLs Jl3|-Jl5[ (which behaves qualitatively as 
well doping), transverse magnetic fields, jlq] and temper- 
ature 0,00 HI in d °p cd SLs - 

Despite its practical and theoretical interest, the ef- 
fect of temperature on electric field domains |2^,^f| and 
current self-oscillations is still poorly understood. Early 
numerical calculations were performed with a fixed drift 
velocity corresponding to a fixed temperature. Us- 
ing the insight provided by these calculations and rea- 
sonable expectations on how drift velocity depends with 



temperature, the fact that oscillatory voltage windows 
widen as the temperature increases was explained. JIJ] 
More detailed experimental studies dealing with the in- 
fluence of temperature on self-oscillations have appeared 
recently. [p7|-^9[ Experimental data show that raising the 
temperature is similar to lowering the SL doping. At low 
temperature a multiplicity of purely static states (corre- 
sponding to coexistence of low and high field domains in 
the SL) was observed. As the temperature increased, 
voltage windows corresponding to self-oscillations ap- 
peared and widened in the SL I-V characteristics. |T|] 
Experimental data were interpreted by using the discrete 
drift model jl3| with a fitted drift velocity. [[19) These au- 
thors concluded that the peak-to-valley ratio in the neg- 
ative differential mobility region of the drift velocity was 
crucial to understand the data. A model including both 
variation of the electron density in the wells and varia- 
tion of the drift velocity with temperature was therefore 
needed. f| 

In a recent paper, we have been able to derive discrete 
drift-diffusion (DDD) models, including boundary condi- 
tions, from microscopic sequential tunneling models, p^ ] 
By using our formulas for the field-dependent drift ve- 
locity at different temperatures (ranging from to 175 
K), we can compare numerical simulations of these sim- 
ple discrete models with Wang et al's experimental data. 
Our results show that increasing temperature facilitates 
current self-oscillations in the second plateau. Further- 
more, our numerical results (based upon microscopically 
calculated drift velocities) agree with the available exper- 
imental data and explain them quantitatively. We ex- 
plain qualitatively why regions of stationary states alter- 
nate with regions of self-oscillations in the temperature- 
voltage phase diagram. Finally and on the basis of our 
numerical simulations, we also explain why the frequency 
may have local maxima in the voltage intervals where 
self-oscillations occur. That the frequency may increase 
with voltage, while the average current simultaneously 
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decreases is thus a consequence of our theory, not an 
anomaly. 0@ 

The rest of this paper is as follows. Section || contains 
a brief description of the DDD model and a calculation of 
its transport coefficients and boundary conditions appro- 
priate for Wang et al's experimental sample. Results of 
numerical simulations of this model and comparison with 
experimental data are reported in Section |ffl| . Section IV 



contains our conclusions. A discussion of the qualitative 
theoretical analysis included in the experimental papers 
is presented in Appendix [A|. 



II. DISCRETE DRIFT-DIFFUSION MODEL 

The main charge transport mechanism in a weakly cou- 
pled SL is sequential resonant tunneling. The character- 
istics of the samples experimentally studied in Ref. |l^] 
are such that the macroscopic time scale of the self- 
sustained oscillations is larger than the tunneling time 
(defined as the time an electron needs to advance from 
one well to the next one) . In turn, this latter time is much 
larger than the intersubband scattering time. Then tun- 
neling across a barrier is a stationary process with well- 
defined Fermi-Dirac distributions at each well. These 
distributions depend on the instantaneous values of the 
electron density and potential drops and vary only on the 
longer macroscopic time scale. The tunneling current 
density across each barrier in the SL may be approxi- 
mately calculated by means of the Transfer Hamiltonian 
method. || The resulting formulas can be used to calcu- 
late the transport coefficients and boundary conditions 
of the following DDD model: @ 

e dF.j n.jv(Fj) , . — it; 



F i -F i _ 1 = -(n i -N%). (2) 

In these equations, e, e and Np are well permittivity, 
minus the electron charge and 2D doping in the wells, 
respectively. C = d + w is the SL period, where d and w 
are the widths of barriers and wells, respectively. Eq. (|l|) 
is Ampere's law establishig that the total current density, 
e J, is sum of displacement and tunneling currents. The 
latter consists of a drift term, eriiv(Fi) / C, and a diffusion 
term, eD(Fi) (nj+i — ni)/ ' C 2 . We have adopted the con- 
vention (usual in this field) that the current density has 
the same direction as the flow of electrons. Eq. ([!]) holds 
for i = 1, . . . , N — 1. Eq. (||) is the Poisson equation, and 
it holds for i = 1, . . . , N. rii is the 2D electron number- 
density at well i, which is singularly concentrated on a 
plane located at the end of the well. Fi is minus an aver- 
age electric field on a SL period comprising the ith well 
and the ith barrier (well i lies between barriers i — 1 and 
i: barriers and N separate the SL from the emitter and 
collector contact regions, respectively). 

Fig. [l] depicts the field-dependent drift velocity at dif- 
ferent temperatures for SL parameter values of Ref. p7| : 



40 periods of 14 nm GaAs and 4 nm AlAs and well dop- 
ing Np = 2 x 10 11 cm~ 2 . It has been calculated from the 
microscopic tunneling current density by the procedure 
explained in Ref. (22). The only adjustable parameter in 
the sequential tunneling formulas is the Lorentzian half- 
width of the scattering amplitudes, 7. To estimate them, 
we have considered that the voltage difference, AV, be- 
tween the peaks of two consecutive branches on the sec- 
ond plateau of the static I—V characteristic is 



AV pa e C 3 - £C2 - 2777. 



(3) 



Here eci is the ith energy level of a given well. For 
7 = 0, the field profile on the second plateau corre- 
sponds to two coexisting electric field domains with fields 
(ec2 - eci)/[e (d + w)] and (e C3 - e C i)/[e (d + w)]. The 
domain walls corresponding to two adjacent branches in 
the I—V diagram are located in adjacent wells. Then the 
voltage difference should be AV w ecz — e C2- In the 
presence of scattering, resonant peaks have finite widths 
which we take as 2777, thereby obtaining (]|). 2t] is an 
adjustable parameter of the order of unity ||. By us- 
ing this formula and the measured current in Ref. [T^ ] 
(Figs. 1, 2, and 3), we find 7 = 18 meV at 1.6 K and 
7 = 23 meV at 140 K for 77 w 0.6. Linear interpolation 
yields the temperature dependence of 7 in the range we 
arc interested in. 
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Electric Field (10 6 V/m) 
FIG. 1. Drift velocity versus electric field for different 
temperatures (starting at K up to 175 K in 25 K steps) 
for a 40- well 14 nm GaAs/4 nm AlAs SL. Well doping is 
N% = 2 x 10 11 cm" 2 . 

Notice that the first peak of the velocity in Fig. |l| 
rapidly disappears as the temperature increases for this 
particular sample. This result might change if we as- 
sume different scattering amplitudes for each of the two 
first subbands of the wells. Moreover, the different ex- 
trema of the velocity curve shift to lower field values as 
the temperature increases. Thus formation of electric 
field domains and current self-oscillations are expected 
for voltages on the second plateau and higher. Multi- 
stable solution branches of the current -voltage character- 
istic curve should also shift to lower voltages and higher 
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currents as the temperature increases, as observed in ex- 
periments. (l9) These effects could not be obtained from 
the fitted drift velocity in Ref. jTjJ. As the diffusion co- 
efficient decreases very rapidly with field, we can safely 
set D = in our DDD model for the experimentally 
observed voltage range. The relevant model is thus the 
well-known discrete drift model of Refs. Jl^Q| with the 
drift velocity deduced from a microscopic calculation of 
the current plus boundary conditions |22| (see Fig. Q). 

To complete the description of our model, we need to 
specify initial, boundary and bias conditions. The dc 
voltage bias condition is 

N 

Cj2 F i = V. (4) 

i=i 

Given that D(Fi) = 0, we need only one boundary con- 
dition specifying Fq (the field at the contact region). We 
will assume that there is an excess electron density in the 
first well due to tunneling from the highly doped contact 
region, § 
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FIG. 2. Bound for the critical doping as a function of 
temperature (full line). Experimental value employed in 
Refs. (rrJ[TJ (dashed line). 



A first interesting conclusion can be drawn effortlessly 
from an analytical upper bound of the critical doping 
above which there are stable static electric field domain 
branches: jll| 



m = (l + c)N 



D ■ 



(5) 



For an appropriately chosen dimensionless positive con- 
stant c, this condition selects recycling of charge 
monopole waves as the mechanism for self-sustained os- 
cillations of the current. The same behavior can be 
obtained from more elaborated boundary condition for 
the microscopic tuneling current between the emitter and 
the neighboring well |22| provided contact doping is suf- 
ficiently high (ohmic behavior). Given the uncertainties 
inherent to contact specification, we have preferred to use 
the phenomenological boundary condition (||) instead. 



III. NUMERICAL SIMULATIONS AND 
COMPARISON WITH EXPERIMENTS 



In this Section, we shall numerically simulate the dis- 
crete drift model for different values of temperature. 
From formula (1) and (2), considering D=0: 



s dF, 
e dt 



v(F) 



Nr 



(Fi 



Ft-i 
N 



)]=J(t), 



i=l 



(0) 
(7) 



The corresponding drift curves are chosen among those 
depicted in Fig. [l] and the boundary condition will be 
(I) with c = 1(T 3 . 



NEc = 



F 



M 



e (dm 



(8) 



In this formula, Fm and F m are the values of the electric 
field which correspond to the maximum and minimum 
of the drift velocity (dm and v m ) on the second plateau. 
The temperature dependence of this critical doping is 
plotted in Fig. ||. We observe that the critical doping 
increases with temperature, indicating that the voltage 
range for which self-oscillations exist increases as tem- 
perature does. In particular, Fig. [2] predicts a transition 
temperature at around 93 K whereas about 140 K is ex- 
perimentally measured. |l7],[l8| Despite the fact that the 
bound (JsJ) is only a rough approximation, the agreement 
is rather good. We shall see below that the agreement 
improves when complete simulations of the DDD model 
are carried out. 

Results of our simulations are presented in Figures |3| 
and U corresponding to a 40-well SL with d = 4 nm, 
w = 14 nm, and N% = 2 x 10 n cm" 2 . Fig. | shows the 
time averaged current-voltage characteristics of such a 
SL for temperatures ranging from 110 to 150 K. In this 
temperature range, there are voltage intervals (with a flat 
form) in which the SL current is stationary interspersed 
in voltage intervals of current self-oscillations correspond- 
ing to recycling of charge monopoles. This agrees with 
experimental results reported by Wang et al. For 
temperatures lower than 110 K or higher than 250 K, the 
ranges of self-oscillations disappear. These figures are 
similar to those reported in Li et al's experiments. |Q9| 
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FIG. 3. I—V characteristics for different temperatures, 
showing stationary (dynamic) states with full (empty) cir- 
cles. The sample is a 40-well 14 nm GaAs/4 nm AlAs SL. 



Well doping is N% 



2 x 10 11 cm" 



10 has been 



used in the numerical simulations. The curve corresponding 
to 150 K has been shifted —0.04 mA for clarity. Lines are 
plotted only for eye-guiding purposes. 

We observe that the I-V curve presents intervals in 
which the average current increases with voltage, fol- 
lowed by intervals in which the average current decreases. 
At lower temperatures the intervals of increasing current 
are wider whereas the opposite occurs at higher tempera- 
tures. Correspondingly, Fig. || shows the frequency of the 
self-oscillations as a function of voltage. The frequency 
of the self-oscillations in such an interval starts increas- 
ing but it drops to a smaller value than the initial one 
at the upper limit of the interval. The amplitude of the 
sclf-oscillations (not shown here) vanishes at the upper 
and lower limits of each voltage interval. This suggests 
that the branches of self-oscillations begin and end at 
supercritical Hopf bifurcations. As the temperature in- 
creases, the region of negative differential mobility in Fig. 
[l] is smoother and the frequency of the self-oscillations in- 
creases (see Fig. |]). In the opposite temperature range, 
at low temperatures, the electric field profiles consist of 
basically two stationary domains joined by a domain wall. 
The I—V characteristic curve has multiple branches cor- 
responding to stationary domains with the domain wall 
located at different wells. This situation resembles that 
obtained as voltage and doping are varied, provided dop- 
ing and reciprocal of temperature are assimilated. In 
Ref. p2| the phase diagram doping-voltage of a doped 
SL was calculated. At low doping (high temperature), 
the electric field inside the SL is almost homogeneous 
and stationary. Above a critical value, branches of self- 
oscillations appear. In this region, there are voltage in- 
tervals of stationary electric field profiles separated by 
intervals of self-oscillations. The latter arise and dis- 
appear (typically) as Hopf bifurcations from stationary 
states. Above a certain doping (low temperature), the 
intervals of self-oscillations vanish and only stationary 
states (consisting of two electric field domains separated 
by a domain wall) remain. Notice in Fig. that there are 



voltage intervals where the oscillation frequency increases 
with voltage, while the average current decreases with 
voltage. This behavior was dubbed anomalous by Wang 
et al, |isf| although it is conveniently explained by the 
discrete drift model equations as shown by our present 
simulations. A qualitative explanation of this behavior 
follows. 
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FIG. 4. Current oscillation frequency versus voltage for 
some dynamic dc bands of the curves shown in Fig. (||). 

First of all, it can be observed in the simulations that 
the maximum current during one oscillation period does 
not vary too much, while the minimum current drops 
precipituously with voltage; see the insets in Fig. ^|. This 
rapid drop of the minimum current and the approximate 
invariance of the maximum current can be understood 
from an asymptotic analysis: || the current approxi- 
mately follows part of the velocity curve, eN^v(F)/C, 
during the motion of a monopole. The monopole is 
created when the current surpasses the maximum value 
eNpV]\i J C. Then the current decreases with time until 
either the monopole exits at the receiving contact or the 
minimum value eNpV rn / C is reached. The first possibil- 
ity is attained at low voltages, the second one at high 
voltages. The current oscillation starts with zero am- 
plitude at the lower end of a voltage interval, so that 
the current does not depart too much from its maximum 
value. As the voltage increases, the minima of the current 
at each oscillation period decrease. What about the de- 
pendence of frequency with voltage? It can be shown Q| 
that the oscillation period T p may be estimated from the 
following formula: 



T p = ——— M — Tf 1 ' 



- 1 



(9) 



Here M is such that N — M represents the number of 
wells traversed by the charge monopole during one oscil- 
lation period, v mon is its average velocity and Tf is the 
monopole formation time (which is an increasing func- 
tion of M [Q). Now, M decreases with voltage owing to 
the dc bias condition, j?]] so that the number of wells tra- 
versed by a monopole grows as the voltage increases (see 
Fig. ||). This means that the first term in ([)]) decreases 
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with voltage and the second one increases with voltage 
(it becomes less negative). Thus there is a competition 
between these two mechanisms: the first tries to make 
the oscillation period decrease with voltage (therefore 
the oscillation frequency increases with voltage), while 
the second term has the opposite effect. Near the low- 
est voltage of an interval for which there are oscillations, 
Vmon ~ vm, and the first term of Eq. ([)]) dominates. As 
the voltage increases, the oscillation amplitude increases 



reported by Wang et al. 17 19| 



and v„ 



decreases, so that the second term becomes 



more important. Of course whether the maximum of the 
frequency is reached immediately or not cannot be said 
from our rough argument. Still, the results of our sim- 
ulations show conclusively that the oscillation frequency 
may reach a single maximum at voltage intervals where 
the current oscillates. 




39 33 
i (index well) 

FIG. 5. Normalized excess of charge at four instants of 
one period of the current self-oscillation at 115 K. Voltages 
are (a) 2.824 V, (b) 2.829 V, (c) 2.837 V and (d) 2.842 V. 
Real-time current traces are plotted in the insets with the 
four instants depicted: ii (full lines), ti (dashed lines), tz 
(dot-dashed lines), and ti (dot-dot-dashed lines). 



IV. CONCLUSIONS 

Starting from a microscopic sequential tunneling model 
of transport in weakly coupled SL, we have derived the 
field-dependent drift velocity for a doped 14/4 SL at tem- 
peratures ranging from to 175K. We found that the first 
plateau rapidly disappears as the temperature increases, 
and therefore we can set the diffusivity to zero in our 
DDD model to study the second plateau. Our numerical 
simulations show that stable solutions change from sta- 
tionary field profiles with two coexisting electric field do- 
mains at low temperature to recycling moving monopoles 
giving rise to current self-oscillations at higher tempera- 
ture. Voltage windows in the I—V diagram appear and 
widen as temperature increases. We observe as well an 
intrincate behavior of the oscillation frequency as a func- 
tion of the DC voltage for different temperatures. These 
findings agree with and explain the experimental data 
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APPENDIX A: WANG ET AL'S THEORETICAL 
EXPLANATIONS 

Theoretical interpretation of experimental results in 
Refs. 0,|l8| is based upon a model proposed earlier 
by X.-R. Wang and Q. Niu (WN). Q Their model is 
mathematically analogous to an earlier model of Laikht- 
man. [E4| It consists of a system of rate equations: 



dt 



k-i{Vi-i) -Ii(Vi), 

Vi - Vi-x = kq u 
N 



(Al) 
(A2) 

(A3) 



with i — 1, . . . N |p5| . Eq. (Al) is charge continuity equa- 
tion for the excess 2D electron charge density at the ith 
well, qi (WN used the notation instead of qi). Ii(Vi) 
is the tunneling current across the ith barrier, which de- 
pends only on the potential drop there, Vi = [ii — (M+ij 
where Hi is the chemical potential at the ith well. (A2) is 
the Poisson equation with k = 4irl 2 /e, I = d + w. Lastly, 



J) is the dc voltage bias condition. The function Ii(Vt) 

is in fact a piecewise linear N-shaped function common 
for all wells, I(Vi); see Fig. 3 of Ref. |§. 

Before commenting on this model, it is convenient to 
simpl ify i t by el imi nating the charge densities, qi. Notice 
that (Al) and (A2) imply the usual Ampere's law 



dt 



h{Vi) 



(A4) 



where i = 0, 1, . . . , N and l(t) is the total cur rent . |2J] 
Thus the model consists of N + 2 equations, (A3) and 
( A4) for N + 2 unknowns, 



Vi, 



0,1, 



, N and 1. The 



current X can be elimina ted by adding all equations ( A4) 
and using that the bias JA3| ) is independent of time. The 
result is the following mean-field model: 



dV 
~dt 



k 



N 



j=0 



v+T g[/i(Vi)-/i(^)]. i = 0,...,N, (A5) 

1 N 

^ftiEw ( A6 ) 



3=0 



5 



Eq. ( A5 ) is reminiscent of the problem of synchronization 
of coupled oscillators of zero frequency at zero tempera- 
ture with the total current playing the role of order, |26[ ] 
parameter. The oscillators Vi try to achieve a stationary 
state such that Ii(Vi) — X, given appropriate proper- 
ties of the functions Ii (positive differential conductivity, 
dh/dVi > 0). H 



1. Critique of the WN model and its analysis 

The main physical objection to the WN model is that 
the sequential tunneling current across a barrier, 7j, de- 
pends explicitly on the electron densities at adjacent 
wells, as well as on electrostatic potentials. This is made 
patent by microscopic derivations, [0,^2) which are con- 
spicuously absent in WN's paper. One unphysical con- 
sequence is that WN's results do not depend explicitly 
on well doping [except that h(Vi) might change in some 
unspecified ad hoc form with doping]. However, experi- 
ments and theory show that stable static electric field do- 
mains are formed at large doping, whereas self-sustained 
current oscillations appear for carrier densities below a 
critical value. |Tl| l Similarly, doping at the injecting con- 
tact (ignored in the WN model) selects the type of charge 
density wave (monopoles or dipoles) responsible for cur- 
rent self-oscillations. Q 

WN's mathematical study of their model contains a 
linear stability analysis of a given unspecified stationary 
state and several unproven statements (some of which 
are even incorrect). Let us be specific. WN claim that 
all eigenvalues A corresponding to inserting Ei = Aie xt 
in the linearized equations are real, which seems reason- 
able. However, later on, they claim that a time periodic 
solution (limit cycle) can appear via a Hopf bifurcation 
from a stationary state. This is an elementary error: a 
necessary condition of a Hopf bifurcation is that a pair of 
complex conjugate eigenvalues cross the imaginary axis. 
The argument they use (illustrated in their Figure 2) to 
"prove" the existence of a limit cycle rests on unproven 
assumptions and, anyway, is not valid in more than two 
dimensions. In fact, they claim that there exists a large 
enough region about an unstable stationary state which 
is invariant under the flow, because the potential differ- 
ence between two adjacent wells cannot exceed the ap- 
plied bias. But no one has proved that this model (not 
to be confounded with physical reality!) posesses such 
desirable property. Furthermore, if there is only one sta- 
tionary state and it becomes unstable by changing a con- 
trol parameter, it will typically do so by having one of 
its eigenvalues changing from negative to positive values. 
(Recall that all eigenvalues are real). Then the bifurcat- 
ing solution will typically (codimension one) be station- 
ary. Thus the situation of WN's figure 2 is unrealistic: 
inside the attractive region depicted, there should be an- 
other stationary (atracting) fixed point. Other more ex- 
otic possibilities are that the flow escapes to decidedly 



unphysical fixed points having some negative Vi or that 
it wanders chaotically between different unstable fixed 
points. 



2. Anomalies 

Wang et al gave an explanation for the fact that 
the self-oscillation frequency may increase with increas- 
ing bias while, at the same time, the mean current de- 
creases. This allegedly anomalous behavior has been ex- 
plained by means of the discrete drift model in Section 
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On the other hand, Wang et al's explanations of 
the "anomaly" are based on the claim that Equations 
(A4) with constant total current X can have a limit cy- 
cle. They also give an estimate of its frequency. These 
arg ume nts are clearly erroneous. In fact, all eq uatio ns 
in (A4) are uncoupled if X is a constant. Then (A4) is 
a one-dimensional autonomous dynamical system, which 
cannot have limit cycles among its solutions. 
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